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EXPERIMENTAL INVESTIGATION OF NOZZLE CRITICAL FLOW CHARACTERISTICS 


OF SODIUM VAPOR 

by Stanley M. Nosek and Louis J. Goldman 
Lewis Research Center 


SUMMARY 

Critical flow rates and pressure ratios of sodium vapor through an axisym- 
metrical convergent -divergent nozzle were measured to determine the character- 
istics of the flow process. The vapor flows were measured from condensate -flows 
by a volumetric measurement technique whereby the condensate, at steady-state 
conditions, was allowed to flow into an auxiliary tank and the rise in level be- 
tween two electrical resistance-type probes was timed. Pressures at the inlet 
and the throat of the nozzle were measured with an inert-gas injection tech- 
nique. 

Flows at nozzle inlet pressures from approximately 5 to 13.5 pounds per 
square inch absolute at saturation temperatures and with up to 100° of superheat 
were investigated. The receiving (condenser) pressure was maintained low enough 
to allow complete expansion through the supersonic nozzle. 

The measured flows without superheat, expressed per unit of nozzle throat 
area and corrected for a flow coefficient C^ of 0.97, could be reasonably ex- 
pressed by w/A-tC^ ss 39.4 p Q where w is the weight flow, A^ is the area of 
the nozzle throat, and p Q is the nozzle inlet pressure. Comparison of the re- 
sults with theoretical one -dimensional flow for an isentropic process obeying 
pv n » c where p is the absolute pressure, v is the specific volume, n is 
the process expansion index, and c is a constant indicated that the saturated 
flows could be represented by an n value of about 1.3. The expansion indexes 
derived from critical pressure ratios agreed favorably with those derived from 
flow measurements. Addition of superheat caused a gradual decrease in flow but 
to the degree applied did not cause a change that might be interpreted as due to 
the point of condensation shifting from upstream to downstream of the throat. 

Comparison of results with theoretical flows for an equilibrium process and 
a supersaturated-frozen process showed closer agreement with those for an equi- 
librium process. A calculation for the process wherein the vapor is supersatu- 
rated but the molecular species are in equilibrium indicated, however, that a 
process of this nature may be the best model for representing the flow. 



IlWRODUCTION 


A knowledge of the flow characteristics of alkali metal vapors under rap- 4 
idly expanding conditions is of practical interest in the design of turbines 
for space power Rankine systems. In particular, it is necessary to have an 
accurate knowledge of the expansion process to properly size the flow passages. 
Without this knowledge, two choices for the process have been usually consid- 
ered, which most probably provide the boundaries of flow predictable by exist- 
ing thermodynamic data. One is the equilibrium expansion process, in which the 
expansion is assumed to occur with equilibrium both between physical phases of 
vapor and liquid and among chemical reaction species of monatomic and poly- 
atomic molecules. The second choice is a supersaturated-frozen expansion in 
which it is assumed that condensation is delayed so that the fluid remains va- 
porous throughout the expansion. It is further assumed that the expansion pro- 
ceeds too rapidly to allow any change in the chemical composition from the ini- 
tial state point (the composition remains frozen) so that the molecular weight 
remains constant. The vapor is presumed to behave as an ideal gas. The theo- 
retically calculated flow between the two processes, as shown in reference 1, 
is significantly different; consequently, a more accurate identification of the 
expansion process is required. 

In references 2 and 3, from flows through nozzles as measured with elec- 
tromagnetic flowmeters, it was indicated that the supersaturated process was 
more plausible for approximating the critical flow of sodium vapor. Both in- 
vestigations, however, left a more accurate identification of the process still 
to be desired. In reference 4, on the basis of accurate pressure measurements 
along a convergent-divergent nozzle and the assumptions that the flow was isen- 
tropic and obeyed the equation pv n = c, it was found that an expansion index 
n of 1.3 best represented the experimental data. Comparison with equilibrium 
and supersaturated processes on the same assumptions would indicate that the 
actual process was between the two. 

To provide further evidence toward identification of the process, the fa- 
cility in references 3 and 4 was modified and operated to obtain more accurate 
nozzle flow measurements. The results are the subject of this report. A vol- 
umetric flow measuring technique was incorporated whereby the condensate was 
allowed, to flow into a tank and the rise between level probes was timed. Flows 
were measured over a range of nozzle inlet pressures from approximately 5 to 13.5 
pounds per square inch absolute at saturation temperatures and with up to 100° of 
indicated superheat. The condenser pressure was maintained from 1 to 3.5 pounds 
■per square inch absolute to provide supersonic flow through the nozzle. 

In addition to the flows, the pressures at the inlet and throat of the 
nozzle were measured with the injection technique described in reference 4. 

The measured flows at saturated conditions are compared with theoretical curves 
of flow calculated on the same assumptions as in reference 4, that is, that the 
vapor flows iscntropieally and obeys pv n = c. Theoretical curves for equilib- 
rium and supersaturated processes were also calculated for comparison. Expan- 
sion indexes were derived from the measured nozzle throat to inlet pressure ra- 
tions and are compared with those determined on the basis of flow. The effect 
of superheat on flow characteristics is also noted. 



SYMBOLS 


A area, sq. in. 

nozzle flow coefficient 

c constant 

g gravitational constant , 32.17 ft/sec 2 

n process expansion index 

p absolute pressure, psia 

T temperature, °F or °R 

v specific volume, lb/cu ft 
w weight flow, lb/hr 

T ratio of specific heat at constant pressure to specific heat at constant 
volume 

Subscripts : 

o nozzle inlet 

t throat of nozzle 


APPARATUS AND IRS TRUMERTAT ION 

The investigation covered by this report was conducted in the sodium flash- 
vaporization facility and research nozzle assembly described in references 3 
and 4. To carry out the research reported herein, the facility was modified by 
adding a condensate flow measuring tank, herein identified as a volumeter. In 
addition, a larger demister was incorporated into the separator, and a swirler 
was inserted into the superheater section between the separator and the nozzle. 


Vapor Generation System 

The vapor generation system is shown schematically in figure 1, where the 
major components and peak design parameters are indicated. The vapor and 
equivalent condensate flow peak value were corrected to show 200 pounds per 
hour rather than 300 pounds per hour as indicated in figure 1 of reference 4 to 
be more consistent with the actual peak value attained. The system design and 
components are described in detail in reference 3. 

The larger demister was deemed to be necessary as a result of the previous 
investigation (ref. 4) in which the dryness of the vapors leaving the separator 
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Nozzle research assembly 



Figure 1. - Schematic of sodium flash-vaporization facility. 



Figure 2. - Swirler and superheater inlet section. 


appeared doubtful. The original demister was a section of 4 -inch pipe filled 
with stainless -steel wire mesh to a depth of approximately 4 inches extending 
into the end cap of the 14 -inch-diameter tank. The end cap was cut off and a 
section of 14-inch pipe, approximately 10 inches in length and completely filled 
with mesh, was added. A new end cap, with the mouth of the 3-inch discharge 
pipe welded flush with the inside surface, was welded on the demister section. 

The swirler was added to augment the demisting and superheat capabilities. 
It was fabricated from 20-gage stainless -steel sheet to form four perpendicular 
vanes and then twisted to provide approximately 270° of swirl. The swirler was 
approximately 2 feet long and fit snugly into the 3-inch discharge pipe from 
the separator. A photograph of the swirler and pipe is shown in figure 2. 
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Research Nozzle 


The nozzle research assembly is included schematically in figure 1. Com- 
plete details of the assembly and pressure measuring technique are presented in 
reference 4. The axial pressure probe was in place , as indicated, but was not 
utilized in this investigation. The nozzle is axisymmetrical and convergent- 
divergent. The convergent section is elliptical in contour and fairs at about 
l/4 inch downstream of the throat into a conical divergent section. At the 
throat, the contour is flat for 0.030 inch. The inlet diameter is 3.06 inches, 
the throat diameter is 0.677, and the outlet diameter is 0.845. The length from 

inlet to outlet is 3.75 inches with the throat 2~ inches from the inlet. The 

area of the throat minus the area of the l/8-inch-diameter probe over the range 
of nozzle inlet temperatures investigated was calculated to be 0.3573 ±0.0006 
square inch. 


Condensate Volumeter 


The condensate volumeter for measuring the vapor flow through the research 
nozzle consisted of a tank nominally 8 inches in diameter and 25 inches long in- 
stalled in parallel with the accumulator as shown in figure 1. Remotely operated 

valves were included, as shown, to manually 
direct the flow as required. A schematic 
drawing of the volumeter with pertinent 
dimensions is shown in figure 3. Installed 
in the volumeter were a low-level and a 
high-level electrical resistance -type probe 
between which the rise in level was timed. 
The volume between the probes as calculated 
from the dimensions, corrected for thermal 
growth at selected condensate temperatures, 
is listed in figure 3. The top of the tank 
was flat to facilitate the installation and 
accurate positioning of the probes. The 
probes were designed on principles set forth 
in reference 5. The electrical resistance 
type was used in preference to the spark 
plug type in order to circumvent the prob- 
lem of shorting out due to condensation on 
the probe. 



Figure 3. - Schematic drawing of condensate volumeter. 


Ins tr ument at i on 

Only instrumentation - temperature, 
pressure, and flow - pertinent to this in- 
vestigation will be presented here. Those 
necessary for operation and control of the 
facility are presented in reference 3. 
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Temperature measurement . - Temperatures were measured with Chromel-Alumel 
thermocouples at locations indicated in figure 1. Two immersed in-well thermo- 
couples were used in the separator, one in the liquid at the bottom, the other 
in the vapor near the top. Both agreed to within 1° or 2° during the course of 
a run. 

The nozzle inlet temperature was measured by a single immersed, sheathed 
thermocouple located approximately 6 inches upstream of the nozzle inlet. The 
absolute value from this thermocouple was found to be in error, and therefore it 
was used only for determining the amount of superheat from changes in its read- 
ings. The separator temperature was used as the nozzle inlet temperature at 
saturated conditions. 

The nozzle discharge temperature was measured in the discharge tee with a 
single immersed, sheated thermocouple. This temperature was measured for com- 
parison with the condenser temperature to note the degree of superheat attained 
by the sudden dissipation of the discharge velocity. In the absence of a qual- 
ity meter, the attainment of superheat at this location would provide some as- 
surance that the quality of the vapors entering the nozzle was high. 

The condenser temperature was measured with two immersed, in-well thermo- 
couples located at the bottom header. The back pressure on the nozzle was con- 
trolled by this temperature, which was always held low enough to assure critical 
flow through the nozzle. The temperature in the volumeter was measured with a 
single, sheathed thermocouple immersed to the depth beyond that of the low-level 
probe as indicated in figure 3. All temperatures were manually recorded as 
readout on a selfbalancing potentiometer indicator calibrated in degrees. 

Pressure measurement . - Static pressures in the nozzle were measured with 
wall taps and the inert-gas injection technique described in reference 4. Nine 
wall taps with 0.020-inch sensing holes were spaced along the nozzle profile. 

In addition, an inlet pressure tap with a 0.060-inch sensing hole was located 
about 6 inches upstream of the nozzle and a discharge pressure tap was located 
in the downstream tee. The throat tap was located at the middle of the 0.030 
flat portion of the profile. Only the inlet pressure and the throat pressure 
were utilized in this investigation. In a few cases, the inlet pressure tap be- 
came plugged so in these cases the pressures between the inlet and the throat 
were used for graphically extrapolating the inlet pressure. The throat tap was 
always clear . 

The pressures were read on a U-tube mercury manometer. The barometric 
pressure was read on-site by evacuating the manometer. The temperature of the 
manometer was also read. 

Vapor flow measurement technique . - The vapor flow through the nozzle was 
measured from its equivalent condensate flow, at steady-state conditions, into 
the volumeter. The value of density that is necessary to convert volumetric 
rate to weight flow rate was calculated from an equation available in reference 
6 or 7. The average temperature of readings taken before, midway, and after a 
run was used. The time for the level to rise between the level probes was meas- 
ured with an electric timer and read to 0.001 of a minute. Starting and stop- 
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ping of the timer was controlled automatically through suitable relays as the 
liquid made contact with the probes. Indicating lights were also automatically 
energized as the liquid reached each level. 

If the accuracies of the volume between the probes, the density, and the 
time measured to fill the volume are considered, the flow rate was believed to 
have been measured to an accuracy well within 1 percent. 


RESULTS OF I3WESTIGATI0N 

The results of the investigation are presented by a compilation of the ex- 
perimental data and a graphical analysis of the results in comparison with theo- 
retical models to determine the characteristics of the flow process through the 
nozzle. With the use of a vapor pressure-temperature relation based on measure- 
ments, the flows at saturated conditions are first correlated against nozzle in- 
let pressure. Comparison of the flow is then made with curves of theoretical 
flow based on assumed processes of equilibrium, supersaturated-frozen, and 
pv n = c. A calculation for a supersaturated-equilibrium process is also in- 
cluded. Expansion indexes n derived for the experimental flows are then com- 
pared with those derived for the experimental critical pressure ratios. The ex- 
pansion indexes from the critical pressure ratios for the equilibrium, 
supersaturated-frozen, and the supersaturated-equilibrium processes are included 
for comparison. Finally, all of the experimental flows are plotted against noz- 
zle inlet temperature to note the effect of superheat. The theoretical trend for 
the bounding equilibrium and supersaturated-frozen processes at a representative 
inlet pressure is used for comparison. 


Experimental Data Obtained 

The experimental data obtained in this investigation are tabulated in 
tables I to III. Table I presents the results of all 102 flow runs, and table II 
presents the averages of the data at each set nozzle inlet condition. Table III 
presents the results of the pressure measurements. The pertinent temperatures 
are included in all tables. The flow rates in tables I and II are expressed per 
unit of throat area and are assumed to be equivalent to vapors entering the noz- 
zle at 100 percent quality. The fact that superheat was attained by the vapors 
at discharge from the nozzle, which may be noted by comparison of the nozzle 
discharge and condenser temperatures, indicated that this assumption is reason- 
able. 


In table III the nozzle inlet pressure, which was measured in inches of 
mercury and corrected for thermal effects, is presented in pounds per square 
inch absolute. The throat pressure is listed indirectly as a ratio to the in- 
let pressure since the primary interest was in arriving at the expansion index, 
which was calculated from 


P 
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TABLE I. - NOZZLE VAPOR FLOW RATE DATA 


Hun 

Sepa- 

rator 

Super- 

heat 

Nozzle 

dis- 

charge 

Condenser 

outlet 

Volu- 

meter 

Vapor 

flow 

rate, 

w/A t , 

lb/(hr ) (sq In. ) 

Run 

Sepa- 

rator 

Super- 

heat 

Nozzle 

dis- 

charge 

Condenser 

outlet 

Volu- 

meter 

Vapor 

flow 

rate, • 
W/A t , 

lb/(hr) (sq In. ) 


Temperature, 1 

°F 




Temperature, T 

, °F 


1 

1433 

0 

1214 

1200 

993 

211.2 

51 

1481 

35 

1248 

1207 

1135 

271.9 

2 

1435 

0 

1210 

1195 

1050 

210.6 

52 

1481 

38 

1250 

1205 

1122 

271.4 

3 

1434 

0 

1210 

1195 

1067 

209.0 

53 

1520 

0 

1268 

1254 

1190 

347.8 

4 

1433 

0 

1207 

1195 

1035 

209.6 

54 

1520 

0 

1270 

1255 

1173 

349.1 

5 

1434 

0 

1210 

1197 

1022 

211.1 

55 

1564 

0 

— 

— 

1230 


6 

1435 

0 

1210 

1197 

1067 

209.2 

56 

1564 

0 



— 

1193 


7 

1434 

0 

1210 

1194 

1017 

216.6 

57 

1565 

0 

1270 

1252 

1237 

431.2 

8 

1435 

0 

1200 

1185 

1085 

209.9 

58 

1565 

0 

1275 

1260 

1250 

424.9 

9 

1520 

0 

1235 

1210 

1161 

339.7 

59 

1565 

0 

1274 

1257 

1245 

428.7 

10 

1520 

0 

1235 

1210 

1153 

340.5 

60 

15 98 

0 

1282 

1252 

1245 

508.0 

11 

1520 

0 

1235 

1211 

1161 

339.5 

61 

1599 

0 

1281 

1252 

1237 

508.9 

12 

1520 

0 

1235 

1210 

1120 

341.0 

62 

1456 

0 

1268 

1259 

1240 

229.0 

13 

1520 

0 

1235 

1210 

1120 

340.7 

63 

1456 

0 

1265 

1259 

1235 

232.0 

14 

1521 

42 

1225 

1204 

1063 

334.1 

64 

1456 

0 

1265 

1255 

1235 

232.4 

15 

1521 

42 

1230 

1205 

1085 

339.7 

65 

1461 

0 

1266 

1255 

1237 

238.7 

16 

1521 

40 

1315 

1302 

1153 

334.7 

66 

1461 

0 

1263 

1253 

1238 

239.7 

17 

1521 

40 

1315 

1300 

1227 

310.8 

67 

1505 

0 

1268 

1255 

1255 

310.7 

18 

1521 

40 

1315 

1301 

1213 

309.1 

68 

1505 

0 

1266 

1254 

1243 

309.8 

19 

1521 

40 

1315 

1301 

1227 

305.6 

69 

1505 

0 

1267 

1254 

1243 

307.6 

20 

1520 

0 

1311 

1300 

1225 

319.3 

70 

1507 

77 

1332 

1273 

1258 

296.2 

21 

1519 

0 

1311 

1302 

1220 

317.8 

71 

1506 

93 

1333 

1267 

1253 

296.5 

22 

1519 

0 

1312 

1302 

1145 

332.5 

72 

1506 

97 

1333 

1265 

1250 

295.8 

23 

1518 

0 

1312 

1300 

1145 

320.0 

73 

1506 

24 

1288 

1265 

1248 

301.4 

24 

1516 

0 

1315 

1304 

1145 

322.4 

74 

1506 

24 

1283 

1265 

1250 

301.4 

25 

1552 

0 

1317 

1305 

1227 

400.3 

75 

1506 

55 

1306 

1264 

1250 

298.2 

26 

1552 

0 

1317 

1306 

1287 

399.3 

76 

1506 

57 

1306 

1264 

1245 

297.6 

27 

1551 

0 

1317 

1305 

1290 

402.0 

77 

1551 

55 

1326 

1269 

1258 

380.3 

28 

1553 

0 

1316 

1304 

1287 

398.5 

78 

1551 

53 

1326 

1269 

1258 

379.9 

29 

1552 

42 

1349 

1302 

1278 

384.1 

79 

1550 

18 

1299 

1270 

1262 

383.5 

30 

1552 

42 

1349 

1305 

1277 

385.6 

80 

1550 

17 

1298 

1269 

1253 

382.0 

31 

1551 

41 

1349 

130!. 

1288 

381.1 

81 

1550 

17 

1297 

1269 

1258 

381.4 

32 

1550 

22 

1333 

1303 

1297 

379.8 

82 

1550 

0 

1284 

1269 

1260 

391.1 

33 

1550 

23 

1335 

1305 

1277 

382.2 

83 

1549 

0 

1284 

1269 

1258 

384.7 

34 

1577 

0 

1320 

1305 

1293 

455.3 

84 

1549 

0 

1284 

1269 

1253 

385.9 

35 

1577 

0 

1320 

1302 

1290 

455.2 

85 

1583 

0 

1280 

1263 

1250 

469.6 

36 

1580 

22 

1337 

1305 

1245 

453. 9 

86 

1583 

0 

1284 

1267 

1252 

469.8 

37 

1580 

22 

1335 

1303 

1240 

454.6 

87 

1461 

0 

1260 

1251 

1240 

238.0 

38 

1578 

0 

1319 

1305 

1245 

464.0 

88 

1460 

0 

1260 

1251 

1240 

238.7 

39 

1578 

0 

1319 

1304 

1245 

462.4 

89 

1460 

0 

1261 

1251 

1238 

236.6 

40 

1600 

0 

1325 

1308 

1245 

515.4 

90 

1459 

23 

1263 

1252 

1240 

229.3 

41 

1600 

0 

1321 

1305 

1245 

518.6 

91 

1459 

23 

1260 

1251 

1238 

227.9 

42 

1600 

0 

1320 

1305 

1245 

515.1 

92 

1460 

23 

1261 

1251 

1240 

227.5 

43 

1565 

0 

1320 

1305 

1300 

424.8 

93 

1461 

60 

1286 

1252 

1238 

227.3 

44 

1565 

0 

1320 

1305 

1280 

424.3 

94 

1460 

61 

1287 

1251 

1238 

226.3 

45 

1565 

21 

1330 

1506 

125 7 

431.0 

95 

1460 

61 

1287 

1251 

1233 

220.7 

46 

1566 

20 

1334 

1305 

1230 

436.3 

96 

1460 

61 

1285 

1250 

1233 

222.3 

47 

1566 

23 

1336 

1305 

1265 

432.1 

97 

1461 

100 

1300 

1252 

1240 

224.3 

48 

1566 

21 

1336 

1306 

1276 

430.8 

98 

1461 

103 

1301 

1252 

1240 

223.3 

49 

1480 

0 

1223 

1206 

1142 

281.5 

99 

1599 

53 

1352 

1286 

1268 

488.7 

50 

1480 

0 

1225 

1207 

1142 

281.8 

100 

1599 

53 

1349 

1282 

1263 

488.4 








101 

1597 

0 

1296 

1281 

1257 

496.6 








102 

1598 

0 

1297 

1282 

1260 

494.7 
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TABLE XI. 


AVERAGE NOZZLE VAPOR PLOW RATE DATA 


Hun 

Sepa- 

rator 

Super- 

heat 

Condenser 

outlet 

Temperature, T, °F 

1 to 8 

1434 

0 

1195 

9 to 13 

1520 

0 

1210 

14 to 15 

1521 

42 

1205 

16 

1521 

40 

1302 

17 to 19 

1521 

40 

1301 

20 to 21 

1520 

0 

1301 

22 

1519 

0 

1302 

23 to 24 

1517 

0 

1302 

25 to 28 

1552 

0 

1305 

29 to 31 

1552 

42 

1304 

32 to 33 

1550 

23 

1304 

34 to 35 

1577 

0 

1304 

36 to 37 

1580 

22 

1304 

38 to 39 

1578 

0 

1305 

40 to 42 

1600 

0 

1306 

43 to 44 

1565 

0 

1305 

45 to 48 

1566 

21 

1306 

49 to 50 

1480 

0 

1207 

51 to 52 

1481 

37 

1206 

53 to 54 

1520 

0 

1255 

57 to 59 

1565 

0 

1256 

60 to 61 

1599 

0 

1252 

62 to 64 

1456 

0 

1258 

65 to 66 

1461 

0 

1254 

67 to 69 

1505 

0 

1254 

70 

1507 


1273 

71 to 72 

1506 


1266 

73 to 74 

1506 


1265 

75 to 76 

1506 


1264 

77 to 81 

1551 


1269 

79 to 81 

1550 

17 

1269 

82 to 84 

1549 

0 

1269 

85 to 86 

1583 

0 

1265 

87 to 89 

1460 

0 

1251 

90 to 92 

1459 

23 

1251 

93 to 96 

1460 

61 

1251 

97 to 98 

1461 


1252 

99 to 100 

1599 


1284 

101 to 102 

1598 


1282 


Measured, 
lb/ (hr) (sq in. ) 

Variation of 
maximum and 
minimum from 
average, 
percent 

Vapor flow rate, w/A t 

210.9 

2.7 


-1.0 

340.3 

.2 


-.2 

336.9 

t .8 

334.7 

0 

308.5 

. 7 


-.9 

318.6 

.2 


-.3 

332.5 

0 

321.2 

+ .4 

400.0 

.5 


-.4 

383.6 

.5 


-.7 

381.0 

±.3 

455.3 

0 

454.3 

±.l 

463.2 

±.2 

516.4 

.4 


-.3 

424.6 

±.l 

432.6 

.8 


-.4 

281.7 

±.l 

271.7 

+ .1 

348.5 

+ .2 

428.3 

.7 


-.8 

508.5 

±.l 

231.1 

.6 


-.9 

239.2 

±.2 

309.4 

.4 


-.6 

296.2 

0 

296.2 

±.l 

301.4 

0 

297.9 

±.l 

380.1 

+ .1 

382.3 

.3 


-.2 

387.2 

1.0 


-.6 

469.7 

0 

237.8 

.4 


-.5 

228.2 

.5 


-.3 

224.2 

1.4 


-1.6 

223.8 

±.2 

488.6 

±.l 

495.7 

±.2 











TABLE III 


NOZZLE PRESSURE DATA 


Sepa- 

rator 

Super- 

heat 

Nozzle 

dis- 

charge 

Condenser 

outlet 

Temperature, T, c 

>F 


inlet throat expansion 

pressure , to inlet index, 
p Q , pressure n^_ 

nsia rat . ic 0 


50A 1480 

50B 1480 

51 1482 

52 1519 

53A 1521 

53B 1520 

59A 1565 

59B 1566 

59C 1566 

61 1598 

66A 1461 

66B 1462 

69 1505 

72 1506 

74 1506 

76 1507 

78A 1550 

78B 1550 

81 1551 

86A 1584 

86B 1585 

89A 1460 

89B 1460 

92 1460 

96A 1461 

96B 1461 

98A 1461 

98B 1461 

100A 1600 

100B 1600 


0 

1225 

0 

1225 

40 

1250 

0 

1268 

0 

1267 

0 

1267 

0 

1270 

0 

1270 

0 

1270 

0 

1279 

0 

1262 

0 

1262 

0 

1266 

97 

1333 

22 

1279 

56 

1308 

56 

1325 

57 

1325 

18 

1296 

0 

1284 

0 

1284 

0 

1260 

0 

1260 

21 

1261 

64 

1287 

64 

1207 

96 

1300 

98 

1300 

54 

1350 

53 

1348 


7.06 
7.02 

a 7 . 01 
a 8 . 82 
8.83 

8.82 

11.18 

11.12 

11.08 

a 13.40 

a 6 . 26 
a 6 . 26 
a 8.07 
a 8.12 
a 8. 02 

8.06 
10.31 
10.30 
10.41 
12.45 

12.38 

6.27 

6.29 

6.30 

6.31 

6.17 

6.20 

6.20 

13.47 

13.53 


0.537 

.539 

.539 

.546 

.548 


1. 

35 

1. 

34 

1. 

34 

1. 

30 

1. 

29 

1. 

30 

1. 

33 

1. 

30 

1. 

29 

1. 

38 

1. 

24 

1. 

23 

1. 

28 

1. 

,32 

1. 

,27 

1. 

,31 

1, 

.33 

1, 

.33 

1, 

.37 

1. 

.36 

1, 

.34 

1 

.27 

1 

.32 

1 

Osl 

ro 

1 

.35 

1 

.30 

1 

.33 

1 

.37 

1 

.41 

1 

.40 


Extrapolated 

















It is assumed that the critical 
pressure of the stream at the wall 
is at the location of the throat 
pressure tap. The critical pres- 
sure ratios are estimated to be 
l/2 to 1 percent accurate, the per- 
cent being better at the higher 
pressures. The expansion index is 
very sensitive to the pressure ra- 
tio, approximately by a factor of 
three, so the accuracy of the ex- 
pansion index is from l|- to 3 per- 
cent. 


Saturated Vapor Pressure - 


Temperature Relation 


The nozzle inlet pressures and 
separator temperatures from table 
III are plotted in figure 4 to ar- 
rive at a calibration curve for use 
in plotting the flow data against 
pressure. This approach was nec- 
essary since the flow and pressure 
could not be measured simultaneous- 
ly and also for the reason that 
both flows and pressures were not measured at all inlet conditions. The results 
also provide a check against vapor pressure curves recommended in references. 
Excellent agreement was found with the equation from reference 8 so that it could 
serve as the calibration curve. In comparison with the curve from reference 7, 
the points fall within 5°, which is within the accuracy, approximately l/2 per- 
cent, of the temperature measurements. In the previous investigation (ref. 4), 
where the nozzle inlet temperature was used, close agreement was found with this 
curve. This provides substantiating evidence that the added demister and swirler 
did not add significant pressure drop to the vapor flow to the nozzle and that 
the separator temperature, therefore, may be used as the saturation temperature 
into the nozzle inlet. 



Weight Flow Correlation 

The average critical flow rates from table II measured at saturated nozzle 
inlet conditions are plotted in figure 5 for correlation against nozzle inlet 
pressure. The pressures were taken from the calibration curve in figure 4. The 
flows are corrected by a flow coefficient to express the results on a theo- 
retical flow basis for subsequent comparisons. The value of was arrived at 

from pretest of the nozzle with air at ambient temperature prior to installation 
in the facility. Over a range of Reynold* s numbers at the throat from 7x10^ to 
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Figure 5. - Variation of experimental critical vapor flow rates with nozzle inlet pressure. 


22X10 4 , the flow coefficient ranged from 0.971 to 0.977. It was estimated, us- 
ing viscosity data from reference 6 for saturated vapor, that the Reynold's num- 
ber would range from 3X10^ to 8X10^ for the conditions investigated herein. A 
constant flow coefficient of 0.97 was therefore considered reasonable to assume. 

It is seen in figure 5 that the least-squares line 


C d\ 


39.4 p 

o 


through the points provides a reasonable representation of the flow. Deviations 
from the line are within ±2 percent at pressures above 9 pounds per square inch 
and up to ±5 percent at pressures below 9 pounds per square inch. The devia- 
tions would have been less, perhaps, had it been possible to factor in the ac- 
tual quality of the vapor. 


Comparison of Experimental Flow with Theoretical Flow 

In figure 6, the flows are compared with theoretically calculated flows. 
Since, at present, there is not a definitive set of thermodynamic data for the 
range of conditions investigated herein, three sources were considered (refs. 7 
to 9 ) for the theoretical calculations. In figure 6(a), properties from ref- 
erence 7 were used. In figure 6 (b), properties from reference 8 were used; and 
in figure 6(c), properties from reference 9 were used. Differences in the sets 
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Figure 6. - Concluded. 

of properties arise from the heat of dimerization value and the vapor pressure 
curve assumed. The heat of dimerization values are -18 200 calories per gram 
mole Nag in reference 7, -16 840 in reference 8, and -17 500 in reference 9. 

The differences in vapor pressure curves are shown in figure 4, where it may be 
noted that the pressures from references 7 and 8 are reasonably close to each 
other and to the pressures measured in this investigation, whereas the pressures 
from reference 9 are substantially lower. 

The theoretical flows shown by the solid lines in figure 6 were calculated 
on the same basis as in reference 4; that is, the flow is assumed to be one- 
dimensional and isentropic and obeys the equation pv n = c. Based on these as- 
sumptions, the following equation for critical flow may be derived: 




300 



2 Nl n +lHn-l) 
rT+I/ 



( 2 ) 


from which a family of curves with selected values of n was calculated and 
plotted as shown. 


Compared with these curves, it is seen in figure 6 that the experimental 
results can be reasonably represented by equation (2) with a constant value of 
n over the pressure range, the value being dependent on the source of proper- 
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ties used. In figure 6(a) , with properties from reference 7, the n value, by 
inspection, is about 1.25; in figure 6(b) with properties from reference 8, it 
is about 1.35; in figure 6(c), with properties from reference 9, it is about 
4 1.33, or nominally 1.3 in all cases. 

Theoretical flows were calculated also for the two processes, described in 
the INTRODUCTION, with are usually considered for tubine design: (l) the equi- 

librium process, and (2) the supersaturated-frozen process. The procedure for 
calculating equilibrium flow rates was the same as that outlined in reference 3, 
that is, finding the maximum ratio of velocity to specific volume along an isen- 
tropic process using tabulated equilibrium thermodynamic properties. The crit- 
ical flow rates for the supersaturated-frozen case were calculated directly 
from equation (2) with the ratio of specific heat at constant pressure to the 
specific heat at constant volume T at the inlet condition used for the expan- 
sion index n. References 7 and 8 do not provide values of IT; therefore, the 
specific heat values in reference 10 for monomer and dimer species were used in 
conjunction with molecular compositions from references 7 and 8 to calculate 
the values. The resulting flows for these two processes are shown in figure 6 
by the dashed lines. In figures 6(a) and (b), the equilibrium process is a 
much closer approximation of the flow; while in figure 6(c), neither one is 
reasonably close. 

Another model that may be considered for comparison is one based on a 
supersaturated-equilibrium process. In this process the vapor expands super- 
saturated with the molecular species in equilibrium so that the molecular weight 
changes during the expansion. A calculation was made in this case with data 
from reference 11 at an inlet temperature of 1625° F (15.18 psia), which is the 
lowest temperature at which properties in this reference are available. An n 
value was derived for this process and the flow was determined from equation 
(2). The result of the calculation is shown in figure 5 (p. 12) where it falls 


, Comparison of Expansion Indexes from 
Flow and Critical Pressure Ratios 


If it is assumed that the least- 
squares line is a better representation 
of the flow than equation ( 1 ) with a 
constant value of n and the least- 
squares line is equated to equation 
(2), values of n as a function of in- 
let pressure may be derived. The re- 
sulting n values were essentially 
linear with inlet pressure and are 
shown by the solid lines in figure 7, 
where the range of values (l.l8 to 
1.28 at 5 psia and 1.3 to 1.4 at 13.5 
psia), again, depends on the properties 


practically on the least-squares line. 


I ! I ! I ! I~T II I I I M l II ! I 

O Derived from critical pressure ratio measurements f 
O Supersaturated-equilibrium process 

Derived from flow rate measurements 

Equilibrium process 

Supersaturated-lrozen process 



6 8 10 12 14 16 

Nozzle inlet pressure, p 0 , psia 

Figure 7. - Variation of expansion index with nozzle inlet pressure. 
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used. The n values calculated from the critical pressure ratio measurements 
are also plotted and substantiate the increasing trend with inlet pressure. 
Agreement is closest with the use of properties from reference 9, although con- 
sidering the accuracy of the data, the agreement is also close to the line • 

formed by using properties from reference 8. The significance of agreement be- 
tween the n values as derived from flow and from critical pressure ratios is 
that it substantiates the reasonableness of assuming that the flow process fol- 
lows pv n = c. The same value of n, thereby, may be used in determining the 
critical flow conditions at the throat that are derivable on this basis (e.g., 
the critical velocity and the critical specific volume) . The use of properties 
from reference 7 necessitates the use of a different value of n for each pa- 
rameter. 

Comparing the results in figure 7 with the expansion indexes for the theo- 
retical equilibrium and supersaturated-frozen processes on the basis of critical 
pressure ratios, shown by the dashed lines, indicates that neither process is a 
reasonable model for the flow. The expansion index for the case of a 
supersaturated-equilibrium process falls reasonably close to the measured values 
and thus provides further evidence that this process may be the best model to 
represent the flow. 


Effect of Superheat 

The purpose of adding superheat in this investigation was primarily to ob- 
serve any pronounced change in the flow that might provide a clue to the nature 
of the process through the nozzle. Based on the analysis in reference 1, ini- 
tially saturated vapors will flow supersaturated to a point in the nozzle where 
condensation will suddenly begin and continue thereafter. The point of conden- 
sation depends on the inlet conditions. Below about 1500° F, the point is 
downstream of the throat and above this temperature is upstream of the throat, 
approaching the inlet as the saturation temperature is increased. Adding super- 
heat at constant pressure reverses the effect, driving the condensation point 
downstream and with sufficient superheat to beyond the throat of the nozzle. 

The effect of the point of condensation having moved from upstream to downstream 
of the throat should be evident by a significant increase in the flow, if it is 
assumed that the process is essentially in equilibrium when the point of conden- 
sation is upstream and is essentially supersaturated-frozen when it is down- 
stream. To determine the magnitude of the change that might be expected, a 
calculation was made at the nozzle inlet pressure corresponding to a saturation 
temperature of 1550° F. At the saturated condition, the supersaturated-frozen 
flow was about 7 percent higher than the equilibrium flow. At 100° of super- 
heat, the flow decreased 1 percent for an equilibrium process and 3 percent for 
a supersaturated-frozen process. The net effect was to reduce the difference 
between the two processes by 2 percentage points so that the difference between 
the two processes at 100° of superheat is, therefore, about 5 percent, which 
should be readily observed. 

The effects of superheat on the two processes are indicated in figure 8 by 
the straight lines originating at the saturated flow curves, which are shown 
for reference. Included in figure 8 are all of the averaged experimental flow 
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Figure 8. - Variation of critical vapor flow with nozzle inlet temperature. 
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data from table II (p. 9). For the experimental points, the nozzle inlet temper- 
ature is the sum of the separator temperature and the amout of superheat meas- 
ured. 


At an inlet condition of 8.8 pounds per square inch absolute, an irregular 
behavior in flow was experienced. With 40° of superheat applied, the condenser 
temperature was raised from 1205° to 1302° F (run 16) . Subsequently, the flow 
decreased with each measurement at fixed inlet condition (runs 16 to 19) . Upon 
removing the superheat, the flow increased but was inconsistent. Later in the 
test program (runs 53 to 54), the flow again showed a wide deviation at this noz- 
zle inlet pressure with no superheat and with the condenser at 1255° F. A con- 
denser temperature between 1302° and 1205° F should not be sufficient to unchoke 
the nozzle, which would explain the marked decrease in flow. Wo unusual behav- 
ior in the facility was observed that might account for the erratic behavior of 
the flow. 

Below and above 8.8 pounds per square inch absolute, where there are suffi- 
cient points, the flow gradually decreases with superheat with no marked increase 
in flow. Either the superheat was insufficient, or the process is other than 
that of the nature assumed, possibly a supersaturated-equilibrium process as was 
shown plausible by the check with saturated flows. 
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SUMMARY OF RESULTS 


Critical flow rates and pressure ratios of sodium vapor flow through an 
axisymmetrical convergent -divergent nozzle were measured to determine the char-* 
acteristics of the expansion process. Flows at inlet pressures from approxi- 
mately 5 to 13.5 pounds per square inch absolute (1425° to 1600° F) and temper- 
atures with up to 100° of superheat were investigated. The measured flows , cor- 
rected for a flow coefficient of 0.97, were compared with theoretical one- 
dimensional isentropic flows obeying pv n = c where p is the absolute pres- 
sure, v is the specific volume, n is the expansion index, and c is a con- 
stant. Comparison of the measure flows with equilibrium and supersaturated- 
frozen flows, as well as a super saturated-equilibrium flow, were also made. Ex- 
pansion indexes derived from critical pressure ratios were compared with the 
values derived from flow. The effect of superheat on the characteristics of 
flow was noted. Three sources of data were considered to investigate the effect 
of differences in available thermodynamic properties on interpretation of re- 
sults. The following results were found: 

1. Over the range of inlet pressures investigated, the critical flows with- 
out superheat and corrected for a flow coefficient of 0.97 could be approximated 
by the straight line w/C^A^ » 39.4 p Q where w is the weight flow, C^ is the 

nozzle flow coefficient, is the area of the nozzle throat, and p 0 is the 

nozzle inlet pressure. Deviations from this line were within ±2 percent at 
pressures above 9 pounds per square inch absolute and as high as ±5 percent at 
pressures below 9 pounds per square inch absolute. 

2. The flows obtained experimentally agreed more closely with theoretical 
flows based on an equilibrium process than on a supersaturated process in the 
range of pressures investigated. A calculation of flow for a supersaturated- 
equilibrirm process agreed very closely with the least-squares line for measured 
flows. The corresponding expansion index was close enough to those from criti- 
cal pressure ratios measured to indicate this process as perhaps the most plau- 
sible. 

3. The expansion index n derived from measured flows and based on 

pv 11 = c increased with nozzle inlet pressure and ranged in values of 1.18 to 
1.28 at 5 pounds per square inch absolute and 1.3 to 1.4 at 13.5 pounds per 
square inch absolute, depending on the source of properties used. A nominal 
value of 1.3 over the range of pressures appeared to be reasonable. The values 
derived from measured critical pressure ratios agreed favorably with those de- 
rived from flow for two of the three property sources and indicated a similar 
increasing trend with inlet pressure. 
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‘ 4. Addition of superheat up to 100° did not cause a change in flow that 
could he attributed to a change in the nature of the process. The flow de- 
creased gradually with increased superheat at any given pressure. 


Lewis Research Center, 

National Aeronautics and Space Administration, 
Cleveland, Ohio, September 13, 1965. 
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